Observation of a New Fluxon Resonant Mechanism in Annular Josephson Tunnel 
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A novel dynamical state has been observed in the dynamics of a perdurbed sine-Gordon system. 
This resonant state, has been experimentally observed as a singularity in the dc current voltage 
characteristic of an annular Josephson tunnel junction, excited in the presence of a magnetic field. 
With this respect, it can be assimilated to self- resonances known as Fiske steps. Differently from 
these, however, we demonstrate, on the basis of numerical simulations, that its detailed dynamics 
involves rotating fluxon pairs, a mechanism associated, so far, to self-resonances known as zero-fleld 
steps. 

PACS numbers: 74.50.-|-r 



Singularities, or resonances, in the current - voltage 
characteristic of a Josephson tunnel junction, reflect the 
underlying dynamics of the phase difference between the 
two superconducting junction electrode order parame- 
ters. Two kinds of singularities are usually observed: 
resonances, known as Fiske steps (FS), that arise when a 
magnetic field is externally applied in the junction plane 
and the so-called zero field steps (ZFS), excited even in 
the absence of an external magnetic field, exclusively in 
large junctions (a Josephson junction is said to be large 
when one of its dimensions is wider than the Joseph- 
son penetration depth \j j]|). Two different resonance 
mechanisms have been proposed to explain the existence 
of these current-voltage singularities: interaction of cav- 
ity modes with the ac Josephson effect and fluxon oscil- 
lations. For small junctions, the theory of Fiske steps 
developed by Kulik [2|, based on the excitation of e.m. 
standing waves, accounts for the experimental observa- 
tions. In large junctions, the fluxon based picture, first 
proposed by Fulton and Dynes 3], constitutes the con- 
vincing explanation frame of reference of all ZFS phe- 
nomenology. In the last approach, fluxons, or particle- 
like magnetic flux quanta (solitons), shuttle back and 
forth along the extended dimension of the junction. Here 
the relevant equation is the one-dimensional sine-Gordon 
equation with appropriate perturbing terms and b.c. [J|. 
There have been attempts to extend the Kulik theory to 
long junctions such that a single type of analysis could 
work for both FS and ZFS. Alternatively, the idea that 
fluxon propagation was responsible also for existence of 
FS, besides ZFS, in long junctions was put forward. This 
hypothesis, sometimes referred as "Samuelsen hypothe- 
sis" la, la] J is based on the simple observation that an ap- 
plied magnetic field renders the junction dynamics asym- 
metric through the boundary conditions, i.e., the fluxon 
propagation becomes unidirectional: fluxons enter one 
of the junction edge and annihilate on the opposite one. 



The situation has been partially redefined by a number 
of experiments performed by Cirillo et al. ^7]. In these 
experiments there is evidence of a possibility which can 
account for both descriptions. There would exist sepa- 
rate regimes in long junctions, depending on the intensity 
of applied magnetic field, in which the two mechanisms 
are well separated and active. At small fields the fluxon 
picture applies, for larger fields, i.e. beyond the threshold 
represented by Hq = 2Xjjc ]^, where jc is the maximum 
critical Josephson current density, the field penetrates 
stably the junction which starts to behave like a short 
junction as far as its magnetic properties are concerned. 
In other words the cavity mode mechanism described by 
Kulik operates beyond Hq while the Samuelsen hypoth- 
esis acts at low field values, when the external magnetic 
field localizes at the edges of the junction. 
Hereafter we report on observations made on an annular 
Josephson junction , i.e., one in which the two electrodes 
are stacked superconducting rings coupled through a thin 
dielectric tunnel barrier. In annular Josephson junctions, 
the above described phenomena have a peculiar charac- 
ter due to the absence of fluxon reflections at the bound- 
aries. Much of the present knownledge about the classical 
fluxon dynamics has been derived by experiments per- 
fomcd on annular Josephson junctions [slllOj. Recently 
Kulik theory has been succefuUy extended to small an- 
nular junctions including 2d effects and the possibility 
of trapped fluxons [llj, |l2|- Moreover quantum effects 
involving trapped fluxons and fluxon-antifluxon (F-AF) 
pair generation have been shown in these devices |l3lll4l |. 
so that fluxon dynamics in annular junctions remains of 
topmost interest in view of their possible use as quantum 
controllable devices ll5|. 

In long annular junctions FS and ZFS appear at the 
asymptotic voltage positions Vn — nc^o/L, where n is 
an integer number, L is the length of the circumference 
of the junction, c is the velocity of light in the June- 



tion (the ultimate fluxon velocity) and $0 the flux quan- 
tum. ZFS dynamics involves free propagation of fluxons 
(antifluxons) around the circle. Unless fluxons are in- 
serted into the junction [1^ |ig, thus changing the so- 
called winding number, ZFS correspond in fact to the 
propagation of fluxon-antifluxon pairs, with zero wind- 
ing number. On the other hand the Samuelsen hyphoth- 
esis for FS in long annular junctions corresponds to the 
following picture: in the presence of an external mag- 
netic field, fluxon-antifluxon pairs are enucleated and suc- 
cessively annihilated in correspondence of two opposite 
points along a diameter normal to the magnetic field di- 
rection where the tangential component of magnetic field 
Ht is maximum and minimum respectively J17l |. These 
two separate mechanisms, in the absence of trapped flux- 
ons, provide through the above formula for Vn the even 
voltage positions (n=2,4,...) for ZFS's and all the posi- 
tions (n=l,2,...) for the case of FS's. 
In this letter we report on experimental and theoreti- 
cal study of phase dynamics underlying resonances ap- 
pearing in the I-V characteristic of a moderately ex- 
tended annular junction. From our analysis a new pic- 
ture of the Samuelsen mechanism appears which permits 
to identify for the flrst time the hybrid character of cer- 
tain fluxon based resonant dynamical states. The un- 
derlying dynamics of the observed states is one in which 
both Samuelsen unidirectional mechanism and ZFS freely 
fluxon propagation coexist. From this point of view they 
stand halfway between FS and ZFS. For larger flelds this 
hybrid dynamics disappears and we are just facing the 
Kulik cavity mode dynamics. 

Experiments have been performed on a Nbhottomi^^^nm) 
/ Albottomi'25nm) /AI2O3/ Altop{2bnm) / Nbtopi^Onm) + 
Nbii;iring{530nm) annular Josephson tunnel junctions. 
Thick Al layers were used in order to reach the low value 
of quality factor, Q, of the junction cavity and, thus, to 
satisfy the low Q condition of the Kulik theory. As for 
the junction geometry, we have chosen island type conflg- 
uration with narrow (« 8 /i m ) wiring leads. The sample 
photograph is shown in the insert of Fig. 1 (a). Internal 
and external junction diameters were 61 ^m and 91 fi m, 
respectively, giving a mean radius oi R ~ 38 fi m. Physi- 
cal parameters of the junction were the following: jc = 54 
A/cm2, Xj = 45 ^nJi, Hq = 0.6 Oe, /" = 27rr = 5.23, 
f = R/Xj ■ The Josephson penetration depth was prac- 
tically equal to the external junction radius so that our 
sample can be considered as a slightly extended annu- 
lar junction. Wc have measured the dependence of the 
amplitudes of resonance steps of the current- voltage char- 
acteristic as a function of parallel magnetic field, H. The 
amplitudes of four resonances were accessible to be reg- 
istred [13. The measurement details are described in 
|l2l |. Fig. 1(a) shows the experimental parallel magnetic 
field dependence of the amplitude of the third resonance 
step (open circles). We used our model theory for small 
annular junctions [12l] as a tool to discriminate between 
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FIG. 1: a) Experimental magnetic field dependence of the cur- 
rent singularity at voltage position corresponding to the third 
resonance step (open circles). The solid line has been obtained 
by an extension of Fiske step Kulik theory to two-dimensional 
annular junction case. Two peaks are visible beside the Kulik 
predicted dependence of the FS. b) Simulated magnetic field 
dependence of the critical current of the third FS obtained 
by numerical solution of Eq.Q. In the simulation I = 5.23, 
Ar = 0.33, a = 1/Q = 0.15. 



features which do not match with the assumption of a 
cavity mode mechanism. The solid curve of Fig. 1(a) is 
the theoretical dependence derived from this theory. The 
experiment shows clearly two peaks, or lobes, in the field 
range 0.46Oe <\ H \< l.QOe which the theory does not 
predict at all. This feature appears in a field region of 
the order of Hq. According to the Samuelsen fluxon pic- 
ture, three pairs should be enucleated at the point where 
Ht is maximum and annihilate on the opposite side. This 
would result in an increment of the phase of 67r in a period 
giving the required FS3 voltage. However the dimensions 
of the junction assume here a critical relevance in that 
the actual realization of the cited mechanism must cope 
with the number of fluxons and antifluxons which actu- 
ally can accomodate into the junction. Differently from 
the case of one-polar states, i.e., only fluxons or only an- 



tifluxons trapped into the junction, existance of stable 
F-AF pairs requires that the junction be sufficiently ex- 
tended to permit the pair dynamics. In particular it ap- 
pears problematic here, for I — 2-Kf < 27r, the occurrence 
of resonances involving more than two pairs [19j . So that 
the above described Samuelsen mechanism requires some 
modifications induced by the junction having such a crit- 
ical dimension. 

In order to investigate deeper this possibility, we 
solved numerically the following 2d perturbed sine- 
Gordon equation which governs the dynamics of the 
phase Lp{r, d, t) in an annular junction: 
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In the above equation lenghts are normalized to the 
Josephson penetration depth Xj and times to the 
Josephson plasma frequency ujj — c/Xj. rie^tji are the 
external and internal magnetic field normalized with 
respect to Xjjc respectively. We assume a 9 dependence 
of normalized bias current which is closer to the exper- 
iment, i.e., two thin leads feed the bias current to the 
junction. In this case (^ = {jb/jc)l: with jf, the bias 
current density and 7 = {■k/A9){1 + ri/2re)Ar, with A6 
the angular width of current lead, Ar the normalized 
junction width, r^ (r^) the normalized external (internal) 
radius. Fig. lb shows the result of simulation of the 
dependence of third resonance current on the external 
magnetic field rje- The internal magnetic field was set to 
0.8577e to optimize the agreement with the experimental 
data, taking into account the screening effects. The 
simulation reproduces fairly well the result of Kulik 
theory for fields higher than Hq. The peaks at small 
fields are also well reproduced. From Fig. lb it is seen 
that these lobes are coexisting with the beginning of the 
subsequent lobes as two separate branches. In fact the 
low field resonant peaks numerical solution branch is 
accessed by FS2 not from the McCumber branch |l| . 
The dynamics involved in the low field peaks is illus- 
trated by the numerical simulation results reported in 
Fig. 2. The phase difference evolution between the 
two electrodes can be followed during a full period 
T = S^o/VsLUj ^ I. The phase distribution appearing in 
Fig. 2 represents only a rough approximation to single 
fluxon propagation, but the overall picture is sufficiently 
clear to draw the following description. A F-AF pair 
(identified with the two points of steepest phase slope) 
is enucleated in the point A and it propagates up the 
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FIG. 2: (Color online) Time plots of the phase at the external 
border of the junction with the corresponding 2d phase snap- 
shots for the annular junctions considered in Fig. l.(^ = 1.104 
and r\t — 1.55 just at the top of the lobe. Eq. (1) is discretized 
in 63x10 spatial sections each long about O.lAj. Section 63 
is followed by section 1. The bias current drive is located just 
at sections 1-5 and 32-36. Point A and B are the maximum 
tangential field locations where pair nucleation (annihilation) 
occurs. Blue arrows indicate the first F-AF pair, red arrows 
the second one. Directions of arrows indicate the fluxon po- 
larity 



positions indicated by the blue arrows as is seen at 
t = 0.000 of Fig. 2. While it moves towards the opposite 
edge (t — 0.664 in Fig. 2 blue arrows) a second pair 
enucleates in the point A and move towards the former 
(see again t — 0.664 in Fig. 2, red arrows). After the 
first pair has reached the point B (t = 1.494 in Fig. 
2) it has an oscillation mode there because energy is 



subtracted and the pair prepares to annihilate. Fig. 2 
(red arrows) at t = 2.998 shows that after the complex 
fourfold coUision {t = 2.158), the second pair survives. 
In fact, it has still enough energy to turn back and 
complete the entire rotation at i = 3.984, there the 
phase flatten and a new pair prepares to enucleate 
at point A (this last is shown at t = 4.980). Things 
can be described as during the collision the second 
incoming pair gains energy at the expenses of the former 
oscillating pair, avoids annihilation, emerges from the 
collision and turns back. From the above analysis we see 
that the mechanism is hybrid: during a single period, 
it involves both a full rotating F-AF pair and a half 
propagating F-AF pair. From this point of view this 
dynamical state shares the nature of a Fiske resonance 
(the second Fiske resonance at V2 in which the two 
pairs both annihilate) with that of a zero field resonance 
(the second zero field resonance at V4 in which the 
two pairs both propagate around the junction making 
a complete turn). This shows that the fiuxon picture 
in the presence of magnetic field is not limited to the 
half propagation mode, on the contrary it can be of 
complex (hybrid) nature. Numerical simulations show 
that hybrid dynamics on third resonance is stable in 
all range of normalized length I between about 4.5 and 
27r. The hybrid lobes for shorter length tend to become 
smaller with respect to Kulik subsequent lobes until they 
disappear when the junction length is less than about 
4.5. On the other hand for I > 2tt the junction dynamics 
switches to Samuelsen mechanism, i.e., three F-AF pairs 
are enucleated on one side and annihilated on the other 
one. In principle hybrid dynamics could exist whenever 
the length of the junction is not sufficient to accomodate 
the required number of pairs to achieve the correct 
voltage. The possibility of hybrid dynamics on other 
resonances, like fourth or fifth, is more questionable and 
maybe depends on the stability of these solutions against 
perturbations induced by other annihilating pairs on the 
ring. 

In conclusion, we have observed a novel fluxon dynamical 
resonant state. In this both half propagating motion of 
a F-AF pair, due to external magnetic field, and propa- 
gation around the whole junction of an other F-AF pair 
coexist. This is probably a further striking observation 
of how nonlinear dynamics described by the sine-Gordon 
equation can be very complex even in a relatively 
simple structure such as an annular Josephson junction. 
It is easy to imagine how hybrid dynamics could be 
generalized to involve both asymmetric propagating and 
back and forth motion of fluxons also in standard long 
linear junctions. Beyond stability the ultimate existence 
of these solutions is tightly related to the dimension of 
their attraction basin which can be difficult to access. 
The possibility of ascribing the observed resonances 
to an hybrid complex dynamics requires a detailed 



analysis of the actual underlaying fluxon propagation 
mode. This should be taken into account specially in 
connection with the potential implementation of vortex 
qubit circuits based on multiply-connected Josephson 
structures [l5j, where a severe control of all possible 
"states" is mandatory. 
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